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Abstract 

The thermodynamic potential of a system of pions and nucleons is computed in- 
cluding the 7rN interactions in the P33 channel. A consistent treatment of the width 
of the resonance in this channel, the A(1232) resonance, is explored in detail. In the 
low-density limit we recover the leading term of the virial expansion for the thermo- 
dynamic potential. An instructive diagrammatic interpretation of the contributions 
to the total baryon number is presented. Furthermore, we examine within a fireball 
model the consequences for the pion spectra in heavy-ion collisions at intermediate 
energies, including the effect of collective flow. A consistent treatment of the A 
width leads to a substantial enhancement of the pion yield at low momenta. 



1 Introduction 



In heavy-ion collisions a zone of hot and dense hadronic matter is created. 
At beam energies on the order of 1 GeV per nucleon or higher, hadron reso- 
nances play an important role in the dynamics of such collisions. Some of these 
resonances have large widths, comparable to or larger than the typical temper- 
ature of the interaction zone. Also stable hadrons generally acquire a width in 
dense matter, due to the interactions with the medium. Consequently, a con- 
sistent theory of the thermodynamics and kinetics of many-body systems with 
short-lived particles is necessarily a basic element in the description of such 
collisions. At present a satisfactory treatment of this problem is still lacking. 

In this letter we explore the consequences of the width of the A(1232) reso- 
nance on the thermodynamics of hot and dense hadronic matter. In heavy-ion 
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collisions at beam energies around 1 GeV per nucleon, the nucleon and the 
pion are together with the A (1232) the most abundant hadrons. By treating 
the A(1232) as an independent field, one effectively accounts for the pion- 
nucleon interactions in the resonant P 33 channel. As a first approximation we 
assume that the interaction zone consists only of these particles. The extension 
of our approach to other channels is straightforward. 

The thermodynamics of such systems has been studied by many authors. In 
calculations based on the Green-function or related techniques the width of 
the A was, so far, either neglected or dealt with in an ad hoc manner [1- 
4]. On the other hand it is well known how to take interactions in general 
and resonances in particular into account within the S-matrix formulation of 
statistical mechanics [5,6]. In section 2 we illustrate the relation between the 
two techniques within a model. 

In order to assess the consequences of the A width for the pion spectrum 
in heavy- ion collisions we consider a simple fireball model (section 3). In this 
model it is assumed that the interaction region, created in such collisions, is in 
local thermal and chemical equilibrium at freeze-out, i.e., when the particles 
decouple. There are two sources for the observed pions: the free pions with a 
thermal distribution present at freeze-out and the pions emanating from the 
decay of thermally excited A-resonances. We show that popular prescriptions 
for evaluating the resonance contribution are incorrect. 



2 Thermodynamic potential and baryon density 



Our starting point for computing the thermodynamic potential Q of a 7rNA-gas 
as a function of temperature T, baryon chemical potential hb and volume V 
is the Green- function approach of Baym and Kadanoff [7]. The temperature 
Green functions are denoted by Qn, Qa for the nucleon and A resonance 
respectively, and by T> n for the pion. The free Green functions are defined by 



g B (p, E m ) = [(iE m - fj, B ) - e B 



T>1 (q, u m ) = (iuj m y - (uj n (q)Y 



e B {p) =m° B +p 2 /2m B (1) 
* (0) = ^Jq 2 + mf , (2) 



where B = N, A while E m and u m denote the Matsubara frequencies for 
fermions and bosons, respectively. We consider an interaction vertex of the 
7riVA type, which connects a pion and a nucleon with a A resonance. The 
following discussion depends only on the topology of the diagrams considered 
and not on a particular choice of the coupling. Hence, we specify the form 
of the interaction only at the end of this section, where numerical results are 
presented. 
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Fig. 1. Diagrammatic representation of (a) the functional and (b) contributions 
to the baryon density. A double line denotes a A, a solid line a nucleon and a 
dashed line a pion Green function. Thick lines correspond to dressed, thin lines to 
free Green functions. The gray circle indicates the baryon number operator. The 
first diagram in (b) yields the nucleon density n^r, the second one the density of 
A resonances in the system while the third one, 5un, corresponds to the 7rN 
component of a physical A (see eqs. (8) and (9)). 



The exact thermodynamic potential £1 of an interacting system of pions, nu- 
cleons and A resonances can be expressed in terms of the full Green func- 
tions [7-9] 



ft = - Tr In 

2 

— Tr In 



- Tr U n V w 



Tr In 



-Qx 1 } - Tr E N G. 



N 



-g^ 1 ] - Tr £ A £ A + $ , 



(3) 



where T, N , X A and 11^ are the self energies of the nucleon, the A resonance 
and the pion. The traces correspond to sums over spin, isospin, and Matsubara 
frequencies as well as integrations over phase space. The functional $ is the 
sum of all two-line irreducible skeleton diagrams where the lines represent 
full Green functions. The self energy of a particle is obtained by a functional 
variation of <3> with respect to the Green function: S B = 5$>/5Qb and TI^ = 
—2 5$/5U n . Diagrammatically, the functional variation is accomplished by 
opening the corresponding line [8]. 

As shown by Baym and Kadanoff [7], the choice of a subset of diagrams for $ 
defines a self-consistent approximation, where every Green function is dressed 
by its self energy obtained by the variation of $. The advantage of such so- 
called conserving or ^-derivable approximations is that they satisfy conser- 
vation laws and guarantee thermodynamic consistency. This approach has in 
the past been applied to a wide variety of many-body problems [9]. For the 
7rNA system the simplest approximation of this type is obtained by choosing 
$ equal to the two- loop diagram shown in fig. 1(a). The corresponding self 
energy of the A resonance is obtained by opening the A line. 



A fully self-consistent treatment of this problem would be quite involved, and 
outside the scope of this exploratory calculation. We therefore relax the self- 
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consistency and consider a low-density expansion^ of eq. (3). In the pion and 
nucleon Green functions we retain only the pole terms, i.e., we use eqs. (1- 
2), albeit with physical masses. On the other hand, for the A, which has a 
large width, the vacuum self energy Sa is treated explicitly. One then finds, 
to leading order in the baryon density 



Q = - Tr In 

2 



Tr In 



-i 



— Tr In 



-1 



(4) 



and QZ 



Q% with m c T[ 



and m° N 



m N , respectively 
Ea- Since the thermodynamic potential in eq. (4) 



where T>\ 

and = (Ql 

is a sum of independent terms involving the vacuum Green functions, it may 
superficially appear that the problem has now been reduced to a trivial one 
with non-interacting pions, nucleons and A resonances. However, as we show 
below, this approximation includes non-trivial interaction effects. In fact, the 
A-resonance term accounts for the 7rN interactions in the P33 channel. 

The traces in eq. (4) can be evaluated using standard techniques (cf. [12]). For 
a spin and isospin symmetric system one finds 



n(T,n B ,V) = 3TV 



-4 TV 



- 16 TV 



d 3 q 
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dE 



B{p,E) In [l + e-^-^l ,(5) 



where the first term is the thermodynamic potential of a free pion gas and the 
second one is that of a free nucleon gas. The last term in eq. (5) accounts for 
the interaction contribution to the thermodynamic potential. It is similar to 
the thermodynamic potential of a free Fermi gas, but involves an integration 
over the energy with a weight function 



B (p, E) = -2 Im 



ggg (P 1 EY ] c r 
dE 



A (p, E) + 2 Im 




(6) 



2 Strictly speaking we loose the nice properties of the conserving approximations 
in this step. However, at least for the problem at hand, the leading term in the low- 
density expansion is close to the self-consistent solution at moderate densities [10] 
(p < 2po). Self consistency is expected to be important, when there is a small scale 
in the problem, see e.g. ref. [11]. 
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Here is the retarded self energy and (G%) = (G° A ) 1 - the retarded 
vacuum Green function of the A resonance, while A denotes its spectral func- 
tion. The functions A and B fulfill the same sum rule, i.e., the integral over 
all energies equals 2n. 

For an interpretation of the weight function B in terms of diagrams it is 
instructive to consider the baryon density 



n B [1 ,/j, B ) =- 



V dfiB 



which can be split into the nucleon contribution 



T,V 



(7) 



n N (T, fi B ) = 4 / ^ ^ N ^\ B) + 1 ( 8 ) 



and the contribution due to the pion-nucleon interactions 

3 

n A (T,u B ) = 16 / / —B(p, E) afJ? 1 , . (9) 

v P ; y (2tt) 3 y 2tt vf ; e^-^ + l 



-co 



Furthermore, as implied by the form (6) for B, the interaction contribution (9) 
to n B consists of two parts. Both of these have a straightforward interpretation: 
The first part, with the spectral function ^4, is the contribution n& of the A 
resonance to the baryon density. The second part, which is proportional to 
dTj^/dE, is the contribution of the nucleons appearing in conjunction with a 
pion in the A self energy. We denote this 7rN pair contribution by 5riN- 

We stress that there is no double counting here. The state vector of an energy 
eigenstate in the interacting system has amplitudes corresponding to a bare 
A and to the non-interacting 7rN states. The spectral function projects the 
interacting states onto the bare A state, while the second term in eq. (6) 
represents the projection onto the 7rN pair states^]. Since the contribution 
of an energy eigenstate to the partition sum depends only on its energy and 
not on the admixture of a particular unperturbed state, both contributions in 
eq. (6) must be included in the thermodynamic potential, eq. (5). 

One advantage of the formulation in terms of Green functions is that dif- 
ferent contributions are easily interpreted in terms of diagrams. The three 

3 A detailed analysis [13] reveals that the second term not only accounts for the 
interacting 7rN pairs, but also subtracts out the contribution of the uncorrelated ones. 
This is necessary to avoid double counting, since the contribution of uncorrelated 
nucleons and pions is taken into account by the first two terms in eq. (5). 
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contributions to the baryon density correspond to the three diagrams shown 
in fig. 1(b). So far the last contribution, 5ri]y, has been neglected in the lit- 
erature on calculations of the thermodynamics of A resonances [2-4]. In such 
calculations an ad hoc prescription was employed, where the spectral function 
A was used in lieu of B. 



In order to make contact with earlier work [5,14], we show that B can be 
expressed in terms of the 7rN-scattering phase shift in the P33 channel. First, 
we note that the direct A-Born graph is the only 7rN-scattering process that 
is generated by the functional fig. 1(a). Hence, the 7rN-scattering amplitude is 
proportional to the vacuum A Green function and the phase shift is given 

byQ 



^33 (E c 



arctan 



Im£g (E c 



— m 



- ReS« (E c 



(10) 



which implies that (cf. eq. (6)) 



B(E cm ) = 2-^5 33 (E cm ) . (li; 



This form of B coincides with that obtained in the S'-matrix formulation of 
statistical mechanics by Dashen, Ma and Bernstein for narrow resonances [5]. 
In the case we consider, namely a broad resonance, there are corrections to this 
expression due to interactions involving more than two particles. However, in 
the low-density limit two-body scattering processes dominate [15]. Thus, when 
the A self energy is properly included, eq. (4) accounts for the 7rN interactions 
in the P33 channel to leading order. When the logarithm in the last term of 
eq. (5) is expanded in powers of exp(— (5{E — /is)), retaining only the lowest 
term, one recovers the well-known expression for the second virial coefficient 
of Beth and Uhlenbeck [14] . 

The expression (11) for the weight function B allows an appealing interpreta- 
tion of eqs. (5) and (9). The energy derivative of the phase shift at the energy 
E cm can be identified with the time delay of a scattered wave due to the 
presence of the scatterer [16]. Hence, the contribution of a resonance to the 
thermodynamic quantities is proportional to the time spent in the resonant 
state [17]. 

The following numerical results were obtained using the 7riVA interaction 
(/tcNa/wO (<f " Sj T a (cf. [18]) with a form factor of the monopole type F(q 2 ) = 

4 Due to Galilean invariance the vacuum Green function and self energy depend 
only on the energy in the irN center-of-momentum frame E cm , i.e. H/\(p,E) = 

cm J ■ 
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Fig. 2. Upper part: Fit of £33 from eq. (10) to the empirical P33 phase shifts [19]; 
Lower part: Weight function B (solid line), see eqs. (6, 11), compared to the spectral 
function A (dashed line). 

A 2 / (q 2 + A 2 ). The parameters are chosen to reproduce the phase shift for 7rN 
scattering in the P33 channel [19] up to a center-of-mass energy of 1.5 GeV 
(upper part of fig. 2). We find f wNA = 3.28, A = 290 MeV and m° A = 1316 
MeV. 

In the lower part of fig. 2 we compare the weight function B with the spectral 
function A. This reveals a qualitative difference close to threshold and a strong 
shift of strength to lower masses, due to the second term of B in eq. (6). The 
weight function B is uniquely determined by the phase shifts, eq. (11), and 
hence a model-independent quantity. On the other hand, the spectral function 
and consequently the separation of %, eq. (9), into contributions from A 
resonances and 7rN pairs are model dependent. 
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3 Pion spectra 



In order to compare our calculation with experimental data we assume that 
at freeze-out the fireball created in a heavy-ion collision is in local thermal 
and chemical equilibrium. Given the population of states ha, which decay 
into a nucleon and a pion after freeze-out, it is straightforward to compute 
the resulting pion spectrum, cf. [3,20]. The pion momentum distribution in 
the A rest frame is folded with the resonance velocity distribution due to 
thermal motion and collective flow. The collective radial flow [4] is taken into 
account by assuming a flow profile, which depends quadratically on the radial 
coordinate^. 

The final 7r° spectrum is then a function of the temperature T, the chemical 
potentials for baryon number and electric charge fiQ as well as the mean 
flow velocity v. The parameters at freeze-out are determined by fitting the ir° 
spectrum [21] of the reaction Au+Au at E^/A = 1 GeV with the constraints 
that the ratio of charge to baryon number equals Nq/Nb = 79/197 and that 
the ratio of tt + mesons to the number of participating baryons is given by [22] 
N + /N B ~ 0.025. 

obs 

In fig. 3 we show the different contributions to the spectrum as a function 
of the transverse pion mass m±_ = \Jp\ + vn\ in a representation, where a 
Boltzmann distribution would appear as a straight line with slope — T~ l . At 
large m±, the spectrum is dominated by thermal pions, while the 7rN pairs give 
an important contribution to the pion spectrum at small transverse masses. 
Since the weigth function B, which determines the sum of the A and 7rN-pair 
contributions, is model independent, the full line is not affected by assumptions 
concerning the structure of the A resonance in vacuum, provided the phase 
shifts are reproduced. At freeze-out we find a temperature of T = 57 MeV, 
a baryon chemical potential of \x& = 830 MeV, a charge chemical potential 
of hq = — 18 MeV and a mean flow velocity of v — 0.36. The corresponding 
freeze-out density is ub = 0.4 p$, where po = 0.16 fm -3 denotes normal nuclear 
matter density. 

As mentioned above the separation of into a bare A and 7rN pair contri- 
bution is model dependent. Nevertheless, it is interesting to see how this is 
realized in a typical case. Using the spectral function A shown in fig. 2 the 
fraction of A resonances at freeze-out is n&/riB = 0.035 while that of A's 
and 7rN pairs is n^/riB = 0.054. When the A resonance is treated as a stable 
particle (Ia = 0) one finds n^/riB = 0.037 for the same values of T and \1b- 



This choice is not crucial, because we find only a weak dependence of the final 
pion spectrum on the assumed analytic form of the flow profile, for a given value of 
the mean flow velocity. 
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Fig. 3. The tt° spectrum of Au+Au at E iab /A = 1 GeV of ref. [21]. The dashed line 
shows the contribution of thermal pions, while the long-dashed line corresponds to 
thermal pions plus those which stem from A resonances, decaying after freeze-out. 
The latter contribution is determined by the spectral function A, fig. 2. Finally, 
the full line includes in addition the contribution of the 7rN pairs. 



In this case the contribution of the correlated 7rN pairs vanishes. 

We note that by arguments analogous to the ones presented here, one finds 
that the invariant mass distribution of 7rN pairs [23] emitted from a fireball 
should not be directly identified with a thermally weighted A resonance. Again 
the correlated 7rN pairs must be accounted for by using B rather than A. 
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4 Summary and conclusions 



We have constructed the thermodynamic potential and the baryon density of 
an interacting 7rN system taking the width of the A resonance consistently 
into account. At low densities, where medium effects can be neglected, the 
weight function for the A resonance is given by its spectral function plus a 
correction term. The latter accounts for the correlated 7rN pair component of 
a physical A resonance. The full weight function can, in agreement with the 
virial expansion, be expressed in a model-independent way in terms of the 
energy derivative of the 7rN scattering phase shift. 

The frequently employed prescription, where the thermodynamics involves a 
trivial spectral function for the nucleon ((^-function) and the free spectral func- 
tion of the A resonance, is incorrect. In such an approach, the contribution 
of 7rN pairs is neglected. These pairs yield an important contribution to the 
population of A-like states and to the low-momentum pion spectrum in rel- 
ativists heavy-ion collisions. When we include all contributions, the thermal 
pions, those due to A decays, as well as the contribution due to 7rN pairs, we 
find a good description of the experimental 7r° spectrum of Au+Au at 1 GeV 
per nucleon. 

We thank E. Grosse, V. Metag, H. Oeschler, P. Senger, C. Sturm and A. Wag- 
ner for rewarding discussions. 
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